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Abstract 

Photon-photon colhsions are investigated in the framework of the two-component Dual 
Parton Model. The model contains contributions from direct, resolved soft and resolved 
hard interactions. All free parameters of the model are determined in fits to hadron- 
hadron and photon-hadron cross section data. The model is shown to agree well to 
hadron production data from hadron-hadron and photon-hadron collisions. The mul- 
tiparticle production in hadron-hadron, photon-hadron and photon-photon collisions as 
predicted by the model is compared. Strong differences are only found as function of the 
transverse momentum variable. The hadron production in photon-photon collisions at 
present and future electron-positron colliders is studied using photon spectra according 
to the equivalent photon approximation, according to beamstrahlung and according to 
backscattered laser radiation. 
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1 Introduction 



The photon, in it's high-energy interactions with hadrons, behaves very much hke a 
hadron, however with cross sections reduced strongly against pure hadronic cross sections. 
Besides this soft hadronic interaction, usually described using the Vector Dominance 
Model (VDM), the photon has a direct (QED) pointlike interaction with the hadronic 
constituents and it has the resolved hard interaction between it's hadronic constituents 
and the hadronic constituents of the target. At moderate energies these hard interac- 
tions of the photons do not change significantly the general picture of photon-hadron 
and photon-photon interactions, this is often forgotten if only the hard part of the pho- 
ton interaction is discussed. Even at high energies, hadronic interactions of photons are 
characterized by soft multiparticle production. Since the soft component of hadron pro- 
duction cannot be understood purely on the basis of perturbative QCD, one has to rely 
on models to calculate the multiparticle final states. The Dual Parton Model (DPM) 
(a recent review is given in Ref. [|^]) has been very successfully describing soft hadronic 
processes in hadron-hadron collisions. Observations like rapidity plateaus and average 
transverse momenta rising with energy, KNO scaling violation, transverse momentum- 
multiplicity correlations and minijets pointed out that soft and hard processes are closely 
related. These properties were understood within the two-component Dual Parton Model 
for hadron-hadron interactions by Aurenche et al. [0, ^ |, |], |], 0, ||. 

Assuming an universal behavior of soft hadronic interactions, it is possible to extent the 
Dual Parton Model to hadronic interactions involving photons. First studies of photon- 
hadron interactions in the framework of the two-component Dual Parton Model were done 
by Engel i, |3> HI- This reaction was studied within other models by various authors, 
one example is the work by Schuler and Sjostrand [|1^, 0. 

Here we apply the model described in |ll]] to the study of hadronic photon-photon 
interactions. In Section II we give a brief introduction to the Dual Parton Model used to 
describe photon-hadron and photon-photon collisions, a complete account of the model 



can be found in An overview of the model realization in the Monte Carlo event 

generator phojet is given. We study with phojet hadron-hadron and photon-hadron 
collisions and compare to data in Section III. In Section IV we compare the properties of 
minimum bias hadron production as calculated with phojet in hadron-hadron, photon- 
hadron and photon-photon collisions. In Section V we discuss the properties of hadron 
production via photon-photon collisions at present and future electron-positron colliders. 
A Summary is given in Section VI. 



2 The event generator PHOJET 

The realization of the Dual Parton Model with a hard and a soft component in PHOJET 
is similar to the event generator dtujet-93 [§, |] simulating p-p and p-p collisions up to 
very high energies. 

In the following, we restrict our discussion of the model to the basic ideas rather than 
giving the complete expressions for all the formulae and quantities entering the model. 



More detailed descriptions of the model are given in H, [11 



In the model, the dual nature of the photon is taken into account by considering the 
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physical photon state as a superposition of a "bare photon" and virtual hadronic states 
having the same quantum numbers J^*" = 1 as the photon. Since the properties of 
the high-mass hadronic fluctuations of the photon are not well known, it is necessary to 
introduce some approximations for calculations. To keep the model as simple as possible, 
two generic hadronic states \qq) and |gg*) have been introduced to describe the hadronic 
piece of the photon. The low- mass state \qq) corresponds to the superposition of the vector 
mesons p, oj and (p and a tt+tt^ background. The state is used as an approximation 
for hadronic states with higher masses (e.g. p', uo' or p"). The physical photon reads 



I7) = |7baro) + |7had) (1) 

with 

g2 g2 e e 

^3 = 1 - 75- - 72- and |7had) = -t~ I^^) + T~ '^^"^ 

Jqq Jqq* Jqq Jqq* 

where e denotes the elementary charge. 

The interactions of the hadronic fluctuations are described within the Dual Parton 
Model in terms of reggeon (M) and pomeron {IP) exchanges. For soft processes, photon- 
hadron duality is used. The energy-dependence of the reggeon and pomeron amplitude 
is assumed to be the same for all hadronic processes. Therefore, data on hadron-hadron 
and photon-hadron cross sections can be used to determine the parameters necessary 
to describe soft photon-photon interactions. However, one does not expect that this 
photon-hadron universality holds for processes involving short distances (high transverse 
momenta). There, long- and short-living fluctuations can contribute. 

In order to use the framework of the QCD-improved parton model with lowest order 



matrix elements [0, |T5[, the pomeron exchange is artificially subdivided into processes 
involving only soft processes and all the other processes with at least one large momentum 
transfer {hard processes). In the model, soft and hard processes are distinguished by 
applying a transverse momentum cutoff to the partons. On Born-graph level, for 

example, the photon-photon cross sections read: 
(i) reggeon and pomeron exchange (soft processes only) 



^Pi') = 9lP,qm 9lP,qq{0) 
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<is) =[ — ] gn,qq{0) 9M,qq{0) ( "j (3) 

with Ajp = ajp(O) — 1 and Ajr = q;jr(0) — 1. Here we denote with a]p{0) (^^(O)) 
the pomeron (reggeon) intercept, and with gr^qq {giR,qq) the couplings of the pomeron 
(reggeon) to the hadronic gg-fiuctuations. 
(ii) hard resolved photon-photon interaction 

QCD 



(4) 
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where f'y,i{xi,Q'^) is the distribution of the parton i in the photon and the sum includes 
all possible parton configurations k, I. 
(iii) single direct interactions 



r da^^^ is i) 

^s-Ms,pT"'') = j dxdiY: fyA^, Q') '"^^ ' ' Q{p^ - pT"^) (5) 

[iv) double direct interactions 



For simplicity, we have written the cross section formulae only for the low-mass state 
\qq). Similar expressions are used for interactions involving the \qq*) state. If not explicitly 
stated, all the calculations have been done using the leading order GRV parton distribution 



functions for the proton and the photon [T^ . 

Assuming Gaussian distributions in impact parameter space, the amplitudes for the 
different processes can be calculated from the cross sections given above. 

The amplitudes corresponding to the one-pomeron exchange between the hadronic 
fiuctuations are unitarized applying a two-channel eikonal formalism similar to |T8|, 0]. In 
impact parameter representation, the eikonalized scattering amplitude for resolved photon 
interactions has the structure 



• / 2 \ 2 

I / e 
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with the eikonal function 

X{s, B) = xs{s, B) + xh{s, B) + B) + xc{s, B). (8) 

Here, Xi{s-,B) denotes the contributions from the different Born graphs: (S) soft part of 
the pomeron and reggeon, (H) hard part of the pomeron (D) triple- and loop-pomeron, 
(C) double-pomeron graphs. To get the photon-photon scattering amplitude, the resolved 
and the direct amplitude given by perturbative QCD are summed up. The complete 
expressions will be given in 

The probabilities to find a photon in one of the generic hadronic states, the coupling 
constants to the reggeon and pomeron, and the effective reggeon and pomeron intercepts 
cannot be determined by basic principles. These quantities are treated as free parameters. 
It was shown in ||^ that it is possible to fix the free parameters by a global fit to proton- 
proton and photon-proton cross sections and elastic slope parameters. In Fig. |l] we show 
the model predictions for the inelastic photon-photon cross section (including quasi-elastic 
vector meson production). The diffractive cross sections of quasi-elastic vector meson 
production (p, uj and 0), single diffraction dissociation and double diffraction dissociation 
are given in Fig. In order to show the strong dependence of the model extrapolations 
on the parton densities of the photon, the cross sections have been calculated with two 
different parametrizations of the parton distribution functions. 
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Once the free parameters are determined, the probabihties for the different final state 
configurations are calculated from the discontinuity of the elastic photon-photon scattering 
amplitude (optical theorem). The total discontinuity can be expressed as a sum of graphs 
with kc soft pomeron cuts, Ic hard pomeron cuts, rric triple- or loop-pomeron cuts, and ric 
double-pomeron cuts by applying the Abramovski-Gribov-Kancheli cutting rules ||19|, pO[] . 
In impact parameter space one gets for the inelastic cross section 

a{kc, Ic, mc, nc, s, B) = — — j j— exp[-2x(s, B)] (9) 

with 

/oo 
d'^B ^(^c, Ic, rric, ric, s, B) ^ dtot - ciqei, (10) 

kc+lc+mc+nc=l 

where cxtot and (Xqei denote the total cross section and the cross section of quasi-elastic 
vector meson production, respectively. We use here the conventions of |^ treating the 
triple- and loop-pomeron cross sections (and hence xd) as negative quantities. In 
the negative sign is explicitly written in the cross section formulae (§). Since the triple- 
, loop-, and double-pomeron graphs are objects involving several pomerons, a further 
resummation is done 0, ^ to allow for the probability interpretation of Eq. (|^). 



In the Monte Carlo realization of the model, the different final state configurations 
are sampled from Eq. (^. For pomeron cuts involving a hard scattering, the complete 
parton kinematics and flavors/colors are sampled according to the Parton Model using 
a method similar to extended to direct processes. For pomeron cuts without hard 
large momentum transfer, the partonic interpretation of the Dual Parton Model is used: 
photons or mesons are split into a quark-antiquark pair whereas baryons are approximated 
by a quark-diquark pair. The longitudinal momentum fractions of the partons are given 



by Regge asymptotics ||22|, E3, |2l, E5|. One obtains for the valence quark (x) and diquark 



'1 — a;) distribution inside the proton 



p{x) ^ - xy-' (11) 

'X 



and for the quark antiquark distribution inside the photon 

p{x) ~ / (12) 
Jx{l — X) 

For multiple interaction events, the sea quark momenta are sampled from a 

p{x) ~ - (13) 

distribution. Note that due to energy-momentum conservation, this distribution is influ- 
enced by the x-distributions of the valence partons and asymmetric multiple interaction 
effects (for example, in photon-proton scattering). This will be discussed in detail else- 
where . The transverse momenta of the soft partons are sampled from the distribution 

--^ ~ exp {-(3p^) . (14) 
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The energy-dependent slope parameter (3 is calculated from the requirement to have a 
smooth transition between the transverse momentum distributions of the soft constituents 
and the hard scattered partons. 

For diffraction dissociation or double-pomeron scattering, the parton configurations 
are generated using the same ideas described above applied to pomeron- 
photon/hadron/pomeron scattering processes. Hence, a diffractive triple-pomeron or 
loop-pomeron cut can involve hard scattering subprocesses resulting in a rapidity gap 
event with jets. According to the kinematics of the triple- or loop-pomeron graphs, the 
mass of the diffractively dissociating systems is sampled from a 1/M^°^^°'' distribution. 
The momentum transfer in diffraction is obtained from an exponential distribution with 
mass-dependent slope (see Ref. 0). For the parton distributions of the pomeron, the 
CKMT parametrization with a hard gluonic component ^ is used. 

Finally, the fragmentation of the sampled partonic final states is done by forming color 
neutral strings between the partons according to the color flow. For soft processes, the 
color flow is approximated using the expansion of QCD for large numbers of colors and 
flavors. This leads to the two-chain configuration characterizing a cut pomeron (as shown 
in Fig. 1^ (a),(b)) and a one-chain system for a cut reggeon. In hard interactions the 
color flow is taken from the matrix elements directly [^]. The leading contributions of 
the matrix elements give a two-chain structure which corresponds to a cut pomeron. For 
example, a cut of a single hard pomeron graph (hard gluon-gluon scattering) is shown in 
Fig. H (c). This method is also applied to the direct photon interactions. 

The chains are fragmented using the Lund fragmentation code jetset 7.3 pO|. 



3 Hadron-hadron and photon-hadron collisions with 

PHOJET 

Hadron production in hadron-hadron collisions has been extensively studied within the 
two-component Dual Parton Model using the dtujet model p, |^ in p-p and p-p colli- 
sions and using the DPMJET-II model in hadron-hadron, hadron-nucleus and nucleus- 



nucleus collisions. Hadron production in photon-hadron collisions is being studied in 
detail using the PHO JET model by Engel 0, |ll| . We present here only a few comparisons 
of PHOJET results with hadron-hadron and photon-hadron data in order to illustrate, 
that the model as formulated in PHOJET is very well able to describe these channels. This 
is certainly required, if we want to apply the model to photon-photon reactions, where 
nearly no experimental data are available. For photon-photon collisions we have to rely 
on the predictive power of the model. 



3.1 Hadron-hadron collisions 

In Figs. I and ^ we compare the transverse momentum distributions and pseudorapidity 
distributions with the results from collider experiments [^, Q . The rise of the plateau 
with the collision energy is understood within the model by one of its most important 
ingredients: the production of multiple soft interactions and multiple hard interactions 
(minijets) rises with energy. The transverse momentum distributions as measured at col- 
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lider energies show besides the soft, exponentially damped, low p± component of hadron 
production the rise of the perturbative hard component with increasing energy. This be- 
comes more and more important with rising energy for minimum-bias hadron production. 

Discussing the changes in hadron production with rising energy, we should indicate, 
that PHOJET in its present version should not be applied for hadron-hadron, photon- 
hadron or photon-photon collisions at collision energies beyond a/s = 1 TeV. One possible 
way, how to apply a model with minijets in the TeV energy region of future proton-proton 
supercolliders has been implemented in dtujet-93 0. This could also be done in phojet, 
however photon-hadron or photon-photon collisions in the multi-TeV energy range seem 
not to be foreseen for the near future. 

In Figs. D and |^ we compare phojet calculations to data on longitudinal distri- 
butions in the fragmentation region using the Feynman scaling variable xp- Feynman 
xjT'-distributions have so far only been measured at fixed target energies, at lab. energies 
well below 1 TeV. It has been shown elsewhere CTI, that the Dual Parton Model shows 



outside the central region (rise of plateau) and the very forward fragmentation region 
(where the diffractive component is dominating) a very good Feynman scaling behavior. 

The seagull-effect, where one plots the average transverse momentum {pj_) as function 
of Feynman xp, demonstrates clearly, that the distributions in transverse momentum 
and longitudinal momentum are correlated in a nontrivial way. There is no factorization 
between the transverse momentum and longitudinal momentum distributions, as it is 
often assumed in oversimplified models of hadron production. In Fig. ^ we compare the 
seagull effect calculated with phojet to data measured inpiab = 360 GeV/c proton-proton 
collisions [133 . 



3.2 Photon-hadron colUsions 



The model for photon-hadron collisions is studied in full detail by Engel Here we 
present only some of this material in order to make the present paper self-contained. 

In Fig. ^ we compare the transverse momentum distribution of charged hadrons cal- 
culated with PHOJET for photon-proton collisions with the photon lab. energy = 140 



GeV with data from the OMEGA Collaboration The data were obtained with a 

tagged photon beam in the energy band 110-170 GeV. The agreement of both distribu- 
tions is excellent, at the largest p± we see already the influence of direct and resolved 
hard collisions. 

In Figs. |l^ and ^ we compare PHOJET results with data from the HERA electron- 
proton collider [R^ . In Fig. ITUI we find a very good agreement with the transverse momen- 



tum distribution of charged hadrons. If we compare to proton-proton collisions (Fig. ^) 
we observe, that the influence of hard collisions in photon-proton collisions is more promi- 
nent than in proton-proton collisions. In Fig. |Tl| we find a good agreement of the inclusive 
charged hadron production cross section as function of the pseudorapidity. We stress, that 
the data as well as the model show a fiat pseudorapidity distribution. 



6 



4 Comparison of minimum-bias hadron production 
in hadron-hadron, photon-hadron and photon-photon 
coUisions 

In this Section we compare the model predictions for inelastic hadron production in 
proton-proton, photon-proton and photon-photon collisions at fixed center-of-mass (CMS) 
energies ^/s. This is the usual way to present data on hadron-hadron collisions at col- 
liders or in fixed target experiments. Collisions of quasi-real photons with protons at 
electron-proton colliders and photon-photon collisions at electron-positron colliders are 
not at fixed photon-proton or photon-photon energy. We will discuss photon-photon col- 
lisions at electron-positron colliders in the next Section. If we consider the three reaction 
channels at fixed energy, then only in order to find the characteristic differences and 
similarities between the three types of collisions. Since elastic hadron-hadron collisions 
usually are excluded studying inclusive secondary distributions, again, in order to find the 
similarities, we also exclude in the reactions with photons the corresponding quasi-elastic 
diffractive channel i.e. 7+7 V + V, (V = p, u, (p) but we include all the other diffractive 
processes. 

We start with the distribution, where the three channels differ most strongly, the 



transverse momentum distribution dN/p±dp±, see Fig. |T2|. The difference is striking. 
The fraction of hard interactions in minimum bias interactions rises from proton-proton 
collisions over photon-proton collisions to photon-photon collisions. The reason for this is 
the direct photon interaction and the fact, that the photon structure function is consider- 
ably harder than the proton structure function. In photon-photon collisions it is easy to 
observe already with moderate statistics hadrons with transverse momenta close to the 
kinematic limit. 

However, these differences in the hard scattering do not strongly influence such average 
properties of the collision as average multiplicities or even average transverse momenta. 
This can be seen from Table 1, where we collect some average quantities characterizing 
nondiffractive proton-proton, proton-photon and photon-photon collisions at CMS ener- 
gies between 10 and 200 GeV. The total and charged multiplicities at all energies are 
rather near to each other in all channels. Probably the differences in the multiplicities 
of newly produced hadrons like vr" and p are more signiflcant, we flnd them at all en- 
ergies rising from p-p over 7-p to 7-7 collisions. Also the average transverse momenta 
rise as expected from p-p over 7-p to 7-7. In Table 1 we give also the number of soft 
pomerons (n)softch. and the number of hard pomerons (n)minijets contributing in average 
to the hadronic flnal state. The numbers given are obtained after kinematical corrections 
due to energy momentum conservation, not the numbers obtained from the unitarization 
step (it turns out, especially at low energies, that for kinematical reasons not all sampled 
pomeron cuts can be generated in the Monte Carlo). At low energy, where the number 
of minijets is very small, we flnd the number of cut soft pomerons to be the same in all 
three channels. The number of cut hard pomerons rises at all energies from p-p over 7-p 
to 7-7. We see also, that the numbers of cut soft and cut hard pomerons are correlated, 
at high energies the number of cut soft pomerons decreases from p-p over 7-p to 7-7. 



In Fig. we compare the longitudinal momentum distributions in the form of XrdN/ dx^ 
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for the three channels. Xr = {p\\/\p\\\)2E/ y/s is the so called radial scaling variable, very 
similar to the Feynman variable xp = 2p\\/y/s. Significant differences between the three 
channels are only found in the region near = 1 or -1. The reason is the single diffrac- 
tive component, which in the p-p and 7-p case leads to the diffractive protons, which are 
obviously missing in 7-7 collisions. 

In Fig. |l^ we compare the seagull effect in the three channels. We find due to the direct 
processes in photon-photon collisions, the rise of {p±) with rising Feynman xp is more 
prominent in photon-photon collisions than in proton-proton collisions. In photon-proton 
interactions we find backwards, in the proton fragmentation region, agreement with the 
p-p collision and in forward direction, the photon fragmentation region, agreement with 
photon-photon collisions. 

In Figs. |1^ and |1^ we plot first logarithmically and then linearly the transverse energy 
distribution dE±/dri. Roughly, these distributions should be equivalent to the pseudo- 
rapidity distribution dN/drj multiplied with the average transverse energy per particle. 
We observe characteristic differences, which can be understood from the features already 
discussed. The transverse energy distribution is wider in photon-photon collisions than 
in proton-proton collisions. The transverse energy distribution rises at all rj from p-p over 
7-p to 7-7. For 7-p the distribution agrees backwards with p-p and forwards with 7-7. 

For p-p or p-p collisions the transverse energy distribution at = is known from ISR 
experiments and experiments at the CERN collider. The values calculated with phojet 
agree well with these measurements. At HERA it was found, that practically the same 
transverse energy at = is found like in p-p interactions and in collisions of real or 



virtual photons with protons This observation agrees well with our results in Figs. 
|15| and 

Finally in Fig. we plot the transverse energy distributions dNjet/dE±jet of jets found 
in the Monte Carlo events from phojet. The jets are searched on hadron level using a 
cone jet finding algorithm with the cone radius R=l. We find again: high ii^_Ljet jets are 
more prominent in 7-7 collisions. 

5 Hadron production in photon-photon collisions at 
present and future electron-positron colliders 

5.1 Photon flux calculation 
5.1.1 Bremsstrahlung 

The fiux of quasi-real photons is calculated using the equivalent photon approximation 
(improved Weizsacker- Williams spectrum |]39| , 0). Within this approximation the 



ep — > eX photoproduction cross section is given by 

da, 



' ep 



dy 

with 



A.ed/) (^7p(^'0) (15) 



27r 



2 ef 1 p2 p2 

min \ min max / 



(16) 
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Here, y and = —p^ denote the energy fraction taken by the photon from the electron 
and the photon virtuahty. aem is the fine structure constant. Taking the kinematic hmit 
-^tnin kin lowcst photou virtuahty allowed one gets with the electron mass me and 



p2 

min,kin 



2 2 
1^ 



(17) 



the spectrum of quasi-real photons 



2ti 



1 + 0^1^(1 



2 

max 



2(1-?/)^ 



y 



ni'ly'^ 



y 



A similar expression holds for ee — ee + X scattering involving quasi-real photons 
only 



dyidy2 



(19) 



5.1.2 Beamstrahlung 

In case of Gaussian beams, the effective beamstrahlung spectrum has been estimated by 
Chen et.al. The dependence of this spectrum on the particle-bunch parameters can 
be expressed by the beamstrahlung parameter Y: 



Y 



(20) 



Here, E denotes the beam energy, Ne is the number of electrons or positrons in a bunch, 
ax and ay are the transverse bunch dimensions, and rg = 2.818 ■ 10~^^mm is the classical 
electron radius. The beamstrahlung spectrum is approximated by ^ 



-1/3 



r(i/3) 

1 — w 



y 



-2/3 (^i _ ^-^y/(i-y) 



+ w 



9{y) 
1 



9iy)Ny 



(1 



with 



~g{y) = l--{l-y) 



2/3 



l-y+{l + y)^l + Y^/^ 



(21) 



(22) 



and K = 2/ (SY), w = 1/ (6y^). The average number of photons emitted per electron 
is given by 

Y 



bal^a^rrie 



2reE Vl + YV3 



(23) 
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5.1.3 Photonemission by laser-backscattering 



Depending on the polarization of the laser light, various photon spectra can be produced 
fi^ , |45| , |46[| . Here we consider only the case of unpolarized laser radiation. Furthermore, 
we assume that the laser frequency is chosen to be below the pair-creation threshold at 
the optimal point. Then, the spectrum of backscattered photons can be written as 

,iaser, ^ _ -0-544 + 2 .17 - 2.63 y + 1.09 



5.2 Hadron production in photon-photon colhsions at electron- 
positron hnear colhders 

There are at present several projects for electron-positron linear colliders under active 
study. Here we pick out for our calculations only one of these projects, the TESLA linear 
collider. Details about the most recent TESLA project were given by Brinkmann . 

Using the formulae discussed in the last subsection we plot in Fig. |l^ the photon spec- 
tra according to the equivalent photon approximation, the beamstrahlung spectrum using 
the bunch parameters as given in the caption of Fig. |18| and a backscattered laser 



spectrum. The photon virtuality was restricted to < 0.01 GeV^/c^. In Table 2 we 
give the average photon-photon energies and the weight factors in /ib for the three photon 



spectra and two energies. From Fig. |I8| and Table 2 we see, that the beamstrahlung spec- 
trum of the TESLA project is the softest of the three photon spectra, the backscattered 
laser spectrum is the hardest. 

Of course, in the case of a linear collider we will always have to consider for background 
problems the superposition of the beamstrahlung spectrum and the equivalent photon 
spectrum. 

A reasonable lowest energy for collisions to be sampled using phojet is a/s^^ = 5 GeV 
(PHOJET would however run without problems even down to ^/s^^ = 2 GeV). Therefore, 
for all applications in this and the following subsection we always cut the photon spectra 
at small y in such a way, that this lower energy cut-off is respected. 

In Figs. |19| and ^ we plot the cross sections adEj_/dri for the transverse energy as 
function of pseudorapidity and da/dr] for the charged hadron production as function of 
pseudorapidity. It is clearly visible, that the backscattered laser spectrum is rather hard 
and has the highest weight. The beamstrahlung spectrum and the equivalent photon 
spectrum are rather comparable, the former has the higher weight, the latter is the harder 
of these two. 

The same differences between the three photon spectra are visible in the cross sections 



da/dp± for charged hadron production as function of the transverse momentum in Fig. ^ 
for the 500 GeV TESLA colhder. 



5.3 Hadron production in photon-photon colhsions at the LEP- 
II electron-positron collider 

We use throughout this section an electron-positron energy of ^/s = 175 GeV for the phase 
II of the LEP collider. Of course here we have only to consider the equivalent photon 
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spectrum. Assuming anti-tagging of the scattered electons and positrons we restrict the 
photon virtuahty to satisfy < 1.75 GeV^/c^. 

The average photon-photon energy at LEP-II with the given electron-positron energy 
will be a/s^^ = 25 GeV. In Table 3 we compare some average properties of the LEP-II 
photon-photon collisions with photon-photon collisions at \/s^^ = 25 GeV. We find, that 
the average properties are rather well represented by photon-photon collisions at ^/s = 
25 GeV. 

In the following we consider a calorimeter detector with the pseudorapidity coverage 
|?7| < 2.1. Applying a jet finding algorithm to this detector we can identify jets in the 
pseudorapidity range |77jct| < l-l- For the jets we use generally a lower jet transverse 
energy cutoff -E±jet > 5 GeV. 

In Fig. 22 we give the jet transverse energy distribution in the acceptance region given 
above together with its decomposition into the different hard interaction mechanisms. 
The single resolved (equal to the single direct) contribution within the acceptance region 
is surprisingly small. At small E±jet the two resolved contribution dominates, at large 
E±]et the direct contribution dominates. 

To study the influence of multiple soft and multiple hard interactions in the model we 
plot in Fig. ^ the jet E± cross section together with further curves, where multiple soft 
or multiple hard contributions or both are rejected. The arbitrary threshold between soft 
and hard partons in the model is p™^°^ = 3 GeV/c. The multiple soft, one hard curve 
coincides nearly with the total E±jet distribution. With the average numbers of minijets 
given in Table 3 the chance to have more than one minijet in one event is small. This is 
different for the multiple soft interactions, with an average (see Table 3) around 1.36 soft 
interactions per event, the chance to have multiple soft chains is rather high. 

To study the influence of an upper visible energy cut we plot in Fig. ^ the jet trans- 
verse energy cross section together with the distribution obtained with the cut -^visible < 
50 GeV within the acceptance region |?7| < 2.1. The 50 GeV cut serves to suppress jets 
coming from decay. For jets with transverse energies below 10 GeV this cut does not 
change the -E±jet distribution drastically. The same is found studying pseudorapidity jet 
profiles with and without the visible energy cut in Fig. Note that the jets become 
more narrow with rising energy and also the pedestal under the jets rises with rising jet 
transverse energy. This again is an effect mainly due to the multiple soft and hard chains. 

Finally, in Fig. ^ we plot the cross section as function of the visible energy within 
the acceptance region. Since we want to search jets with a -E±jct > 5 GeV cutoff, we 
use a lower threshold for the photon-photon energy ^/s^^ = 10 GeV. The plot gives the 
visible energy distribution without restrictions and in a second curve the visible energy 
obtained if we demand at least one jet with E^j^t > 5 GeV and impose the visible energy 
cut -^visible < 50 GeV within the acceptance region given above. We observe: nearly each 
event with a visible energy above 30 GeV contains at least one jet with E±jct ^ 5 GeV. 



6 Conclusions and summary 

The PHOJET model can be used to calculate hadronic events in hadron-hadron, photon- 
hadron and photon-photon collisions. The model is found to agree well with data in 
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hadron-hadron and photon-hadron collisions, the predictions for photon-photon collisions 
do not need any new parameters. 

Multiple soft and multiple hard interactions (minijets) lead to a rise of the rapidity 
plateau, which agrees in hadron-hadron and photon-hadron collisions very well with the 
rise of the plateau observed experimentally. 

Minimum bias hadron production in hadron-hadron, photon-hadron and photon-photon 
collisions of the same CMS energy is remarkably similar. To see this, one has to restrict 
the comparison to inelastic events and to exclude also the diffractively produced vector 
mesons in reactions involving photons. The only striking differences appear in the trans- 
verse momentum distribution or distributions, where the transverse momentum behavior 
is essential. This difference can be understood to be due to the direct photon interaction 
contribution and due to the photon structure function being considerably harder than 
hadronic structure functions. 

PHOJET can be applied also to photon-photon interactions in electron-hadron and 
electron-electron or electron-positron colliders. Photon spectra according to the equiva- 
lent photon approximation, according to beamstrahlung spectra and backscattered laser 
radiation in linear electron-positron colliders are implemented at present. It is easy to 
implement other similar photon spectra. With these possibilities PHOJET should be a 
tool suitable to study the photon-photon background to other interesting reactions at 
such colliders. 
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Table 1 Comparison of average quantities characterizing hadron production in non- 
diffr active p-p, 7-p and 7-7 collisions at cms energies between 10 and 200 GeV. All 
average transverse momenta are given in GeV/c. 





Quantity 


P-P 


7-p 


7^ 


10 




11.2 


11.1 


11.7 


10 


rich 


6.65 


6.53 


6.86 


10 




2.17 


2.44 


2.88 


10 


"p 


0.027 


0.063 


0.11 


10 


^ Pi-ch ^ centr.rj 


0.39 


0.38 


0.42 


10 


< P± >7r- 


0.32 


0.33 


0.36 


10 


< P± >p 


0.41 


0.43 


0.47 


10 


< fl >softch. 


1.16 


1.19 


1.22 


10 


< n ^ mini jets 


0. 


0.00004 


0.0020 


20 


ritot 


16.4 


16.6 


17.1 


20 


rich 


9.64 


9.71 


10.00 


20 


n-K- 


3.44 


3.78 


4.18 


20 


rip 


0.086 


0.14 


0.20 


20 


^ PJ-cH ^centr.T] 


0.37 


0.38 


0.44 


20 


< P± >7r- 


0.32 


0.34 


0.38 


20 


< P± >p 


0.42 


0.45 


0.52 


20 


< ri >softch. 


1.26 


1.33 


1.31 


20 


^ n ^ mini jets 


0.0003 


0.0025 


0.028 


50 


ritot 


24.8 


26.5 


26.9 


50 


rich 


14.5 


15.5 


15.6 


50 


n-K- 


5.49 


6.19 


6.53 


50 


rip 


0.21 


0.27 


0.34 



16 



Continuation of Table 1 
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0.38 


0.42 


200 


< P± >p 


0.47 


0.53 


0.64 


200 


< ^ >softch. 


1.59 


1.87 


1.29 


200 


^ ^ mini jets 


0.17 


0.36 


1.01 
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Table 2 Average photon-photon energies and weight factors in /i6 for the three photon 
spectra and two hnear coUider energies. 





Photon spectrum 




Weight (lib) 


500 


W.-W. 


50 


0.0067 


500 


Beamst. 


16.7 


0.038 


500 


B. Laser 


252 


0.47 


1000 


W.-W. 


105 


0.0075 


1000 


Beamst. 


41 


0.026 


1000 


B. Laser 


509 


0.59 



Table 3 Comparison of average quantities for LEP-II photon-photon coUisions with 
photon-photon coUisions at \/s = 25 GeV. 



Quantity 


LEP-II 77 


77 at 25 GeV 


ntot 


15.6 


17.2 


rich 


9.19 


10.1 


"'7r+ 


3.74 


4.12 




3.74 


4.12 


np 


0.17 


0.19 


< P± >7r- 


0.37 


0.37 


< P± >p 


0.50 


0.50 


< n >softch. 


1.36 


1.44 


^ n ^ mini jets 


0.035 


0.028 


Weight (fjib) 


0.0032 
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Figure 1: Inelastic photon-photon cross sections as calculated in the model compared 
with experimental data at low energies E3, E^, [ 



51, 521. The two curves from the 



model were calculated using the GRV LO photon structure function and the SaS 
2M photon structure function [^. The differences between both curves at high energy 



demonstrate the uncertainties of the predictions due to the limited data available on the 
photon structure function. Our curve calculated with the SaS 2M structure function 
agrees practically with the cross section calculated with the same structure function but 
using another model by Schuler and Sjostrand |^ . 



Figure 2: Diffractive cross sections as calculated with PHOJET using the GRV LO photon 
structure function and the SaS 2M photon structure function |BB[ . The upper curve is 



for each of the three cross sections the one obtained with the GRV LO structure function. 



Figure 3: Unitarity cut of a one-pomeron graph: the unitarity sum including all possible 
final states is subdivided into final states with low-p_L partons and into final states with 
at least one parton satisfying > p™^°^. 



Figure 4: Comparison of transverse momentum distributions of charged hadrons with 
collider data at ^/s = 200 GeV The calculation uses the Dual Parton Model code 

PHOJET. 



Figure 5: Pseudorapidity distributions of charged hadrons produced in p-p collisions 
as calculated with phojet are compared to collider data from the UA-5 Collaboration 
3|, g for the energy = 200 GeV. 



Figure 6: Comparison of Feynman xp distributions of 7r+-mesons produced in proton- 
proton collisions at 205 and 175 GeV. The experimental data are from Kafka et al. 
and from Brenner et al. The data from both experiments agree rather well with 

each other. The calculation uses the Dual Parton Model phojet. 



Figure 7: Comparison of Feynman xp distributions of vr'-mesons produced in proton- 
proton collisions at 400 and 360 GeV. The experimental data are from the Aguilar-Benitez 
57[] and from the EHS-RCBC Collaboration The data from both experiments 



et al. 

agree rather well with each other, in fact most of the data points of [ESI are below the 



57[] data. The calculation uses the Dual Parton Model phojet. 
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Figure 8: Comparison of the seagull effect in the reaction p+p— h +X at 360 GeV. The 



data are from the EHS-RCBC Collaboration The calculation uses the Dual Parton 
Model PHOJET. 



Figure 9: Comparison of the transverse momentum distribution of charged hadrons with 
< Xi^' < 1 as calculated with phojet at the average photon energy of 140 GeV with 
data. The data are from the OMEGA Collaboration ||36[ measured with a tagged photon 
beam in the energy band 110-170 GeV. 



Figure 10: Comparison of the transverse momentum distribution of charged hadrons 
produced by quasi-real photons in electron-proton collisions at HERA [pTf with the pho- 



jet calculation [jTO[. The data and the calculation is in the pseudorapidity region -1.5 

< '7lab < 1-5. 



Figure 11: Comparison of the inclusive pseudorapidity cross section of charged hadrons 



produced by quasi-real photons in electron-proton collisions at HERA with the pho- 
jet calculation |l^ . 



Figure 12: We compare at the collision energy ^/s = 200 GeV the transverse momentum 
distribution in invariant form for all charged hadrons produced in proton-proton, photon- 
proton and photon-photon collisions. The calculation was done with phojet for inelastic 
collisions, excluding in photon-proton collisions the V-p and in photon-photon collisions 
the V-V diffractive production of vector mesons V = p, u and 0. 



Figure 13: We compare at the collision energy ^/s = 200 GeV the radial Xr = lEj^ 
distribution in invariant form for all charged hadrons produced in proton-proton, photon- 
proton and photon-photon collisions. The calculation was done with phojet for inelastic 
collisions, excluding in photon-proton collisions the V-p and in photon-photon collisions 
the V-V diffractive contributions. 



Figure 14: We compare at the collision energy ^/s = 200 GeV the average transverse 
momentum of charged hadrons produced in proton-proton, photon-proton and photon- 
photon collisions as function of the Feynman xp variable (seagull effect). The calculation 
was done with phojet for inelastic collisions, excluding in photon-proton collisions the 
V-p and in photon-photon collisions the V-V diffractive contributions. 
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Figure 15: We compare at the collision energies ^/s = 20 GeV and 200 GeV in a log- 
plot the distribution of the transverse energy Ej_ as function of the pseudorapidity t] for 
proton-proton, photon-proton and photon-photon collisions. The calculation was done 
with PHOJET for inelastic collisions, excluding in photon-proton collisions the V-p and in 
photon-photon collisions the V-V difFractive contributions. 



Figure 16: We compare at the collision energies y/s = 20 GeV , 200 GeV and 1000 GeV in 
a linear plot the distribution of the transverse energy E± as function of the pseudorapidity 
r] for proton-proton, photon-proton and photon-photon collisions. The calculation was 
done with phojet for inelastic collisions, excluding in photon-proton collisions the V-p 
and in photon-photon collisions the V-V diffractive contributions. 



Figure 17: We compare at the collision energy ^/s = 200 GeV the transverse energy dis- 
tribution for hadronic jets (identified using a jet-finding algorithm) produced in proton- 
proton, photon-proton and photon-photon collisions. The calculation was done with PHO- 
JET. 



Figure 18: Photon fluxes at a ^/s = 500 GeV linear collider TESLA [^]. Given are the 
improved Weizsacker- Williams spectrum, the beamstrahlung spectrum using the bunch 
parameters Ne = 1.8 10^°, 
backscattered laser spectrum. 



598 nm, cr,. 



6.5 nm and cr? = 0.5 mm [H7|, and a 



Figure 19: 

in fib GeV at the y/s 



Cross section weighted transverse energy distributions adE^/drj measured 

p. 



Given are the 



500 and 1000 GeV TESLA linear colliders 
distributions for the Weizsacker- Williams photon spectrum, the beamstrahlung spectrum 
and a backscattered laser spectrum. 



Figure 20: Pseudorapidity cross sections da/dr], measured in /ib per pseudorapidity 
unit, at the y/s = 500 and 1000 GeV TESLA linear colliders g^l- Given are the distri- 
butions for the Weizsacker- Williams photon spectrum, the beamstrahlung spectrum and 
a backscattered laser spectrum. 



Figure 21: Transverse momentum cross sections da/dp± at the y/s = 500 GeV TESLA 
linear collider [^]. Given are the distributions for the Weizsacker- Williams photon spec- 
trum, the beamstrahlung spectrum and a backscattered laser spectrum. Please note that 
the p± distributions for the Weizsacker- Williams spectrum and the beamstrahlung spec- 
trum cross. At low p± the beamstrahlung dominates, at high p± the W.-W. spectrum 
dominates. 



Figure 22: Decomposition of the jet cross section dcrjet/rf-Exjet as function of the jet 
transverse energy. The jets with -E±jet > 5 GeV were found from the phojet events using 
a cone algorithm. The jet acceptance region was restricted to Ir^jctl < 1.1- 
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Figure 23: Jet cross section da^et/ dE±^f.t as function of the jet transverse energy. In order 
to study the influence of muhiple soft and hard interactions in the model, we consider 
restricted cases with only single soft or single hard collisions. The curve with multiple 
soft, single hard collision coincides nearly with the full curve. The jets with -E_Ljet > 5 
GeV were found from the phojet events using a cone algorithm. The jet acceptance 
region was restricted to |?7jet| < 1.1. 



Figure 24: Jet cross section dajet/dE±jet as function of the transverse energy. Here we 
give in addition the cross section with a visible energy cut (inside a calorimeter with \ri\ < 
2.1 the cut is -Evisibie < 50 GeV). The jets with E±jet > 5 GeV were found from the phojet 
events using a cone algorithm. The jet acceptance region was restricted to |?7jct| < l-l- 



Figure 25: Jet profiles in the pseudorapidity variable with and without the visible energy 
cut described in the caption of Fig. 



Figure 26: Cross section as function of the visible energy inside a calorimeter with 
\ri\ < 2.1. Only events with y/s^^ > 10 GeV are sampled from the improved Weizsacker- 
Williams spectrum. In a second curve we request jets with -E_Ljct > 5 GeV and impose a 
visible energy cut of 50 GeV. 
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Figure 1: Inelastic photon-photon cross sections as calculated in the model compared 
with experimental data at low energies [48, 49, 50, 51, 52]. The two curves from the 
model were calculated using the GRV LO photon structure function [17] and the SaS 
2M photon structure function [53]. The differences between both curves at high energy 
demonstrate the uncertainties of the predictions due to the limited data available on the 
photon structure function. Our curve calculated with the SaS 2M structure function 
agrees practically with the cross section calculated with the same structure function but 
using another model by Schuler and Sjostrand [54]. 
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Figure 2: Diffractive cross sections as calculated with PHOJET using the GRV LO photon 
structure function [17] and the SaS 2M photon structure function [53]. The upper curve is 
for each of the three cross sections the one obtained with the GRV LO structure function. 




Figure 3: Unitarity cut of a one-pomeron graph: the unitarity sum including all possible 
final states is subdivided into final states with low-p± partons and into final states with 
at least one parton satisfying pj_ > p'^^^°^. 
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Figure 4: Comparison of transverse momentum distributions of charged hadrons with 
coUider data at -^i = 200 GeV [32]. The calculation uses the Dual Parton Model code 
PHOJET. 
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Figure 5: Pseudorapidity distributions of charged hadrons produced in p-p collisions as 
calculated with PHOJET are compared to collider data from the UA-5 Collaboration [33, 
34] for the energy ^/s = 200 GeV. 
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Figure 6: Comparison of Feynman xp distributions of tt "'"-mesons produced in proton- 
proton collisions at 205 and 175 GeV. The experimental data are from Kafka et al. [55] 
and from Brenner et al. [56]. The data from both experiments agree rather well with 
each other. The calculation uses the Dual Parton Model PHOJET. 
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Figure 7: Comparison of Feynman xp distributions of 7r~-mesons produced in proton- 
proton collisions at 400 and 360 GeV. The experimental data are from the Aguilar-Benitez 
et al. [57] and from the FHS-RCBC Collaboration [58]. The data from both experiments 
agree rather well with each other, in fact most of the data points of [58] are below the 
[57] data. The calculation uses the Dual Parton Model PHOJET. 
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Figure 8: Comparison of the seagull effect in the reaction p+p— > h~ +X at 360 GeV. The 
data are from the FHS-RCBC Collaboration [35]. The calculation uses the Dual Parton 
Model PHOJET. 
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Eigure 9: Comparison of the transverse momentum distribution of charged hadrons with 
< a^F ^ 1 as calculated with PHOJET at the average photon energy of 140 GeV with 
data. The data are from the OMEGA Collaboration [36] measured with a tagged photon 
beam in the energy band 110-170 GeV. 
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Figure 10: Comparison of the transverse momentum distribution of charged hadrons pro- 
duced by quasi-real photons in electron-proton collisions at HERA [37] with the PHOJET 
calculation [10]. The data and the calculation is in the pseudorapidity region -1.5 < ?7 lab ^ 
1.5. 
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Figure 11: Comparison of the inclusive pseudorapidity cross section of charged hadrons 
produced by quasi-real photons in electron-proton coUisions at HERA [37] with the PHO- 
JET calculation [10]. 
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Figure 12: We compare at the collision energy = 200 GeV the transverse momentum 
distribution in invariant form for all charged hadrons produced in proton-proton, photon- 
proton and photon-photon coUisions. The calculation was done with PHOJET for inelastic 
collisions, excluding in photon-proton collisions the V-p and in photon-photon collisions 
the V-V diffractive production of vector mesons V = co and (p. 
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Figure 13: We compare at the collision energy = 200 GeV the radial = 2E/y^ 
distribution in invariant form for all charged hadrons produced in proton-proton, photon- 
proton and photon-photon collisions. The calculation was done with PHOJET for inelastic 
collisions, excluding in photon-proton collisions the V-p and in photon-photon collisions 
the V-V diffractive contributions. 
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Figure 14: We compare at the collision energy = 200 GeV the average transverse 
momentum of charged hadrons produced in proton-proton, photon-proton and photon- 
photon collisions as function of the Feynman xp variable (seagull effect). The calculation 
was done with PHOJET for inelastic collisions, excluding in photon-proton collisions the 
V-p and in photon-photon collisions the V-V diffractive contributions. 
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Figure 15: We compare at the collision energies = 20 GeV and 200 GeV in a log- 
plot the distribution of the transverse energy E± as function of the pseudorapidity rj for 
proton-proton, photon-proton and photon-photon collisions. The calculation was done 
with PHOJET for inelastic collisions, excluding in photon-proton collisions the V-p and in 
photon-photon collisions the V-V diffractive contributions. 
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Figure 16: We compare at the collision energies ^/s = 20 GeV , 200 GeV and 1000 GeV in a 
linear plot the distribution of the transverse energy E± as function of the pseudorapidity 
T] for proton-proton, photon-proton and photon-photon collisions. The calculation was 
done with PHOJET for inelastic collisions, excluding in photon-proton collisions the V-p 
and in photon-photon collisions the V-V diffractive contributions. 
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Figure 17: We compare at the collision energy -^s = 200 GeV the transverse energy dis- 
tribution for hadronic jets (identified using a jet-finding algorithm) produced in proton- 
proton, photon-proton and photon-photon collisions. The calculation was done with PHO- 
JET. 
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Figure 18: Photon fluxes at a ^/s = 500 GeV linear collider TESLA [47]. Given are the 
improved Weizsacker- Williams spectrum, the beamstrahlung spectrum using the bunch 
parameters = 1.8 10^°, (t^^ = 598 nm, (Ty = 6.5 nm and cr^ = 0.5 mm [47], and a 
backscattered laser spectrum. 
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Figure 19: Cross section weighted transverse energy distributions crdE±/dr] measured in 
/ib GeV at the = 500 and 1000 GeV TESLA hnear coUiders [47]. Given are the 
distributions for the Weizsacker-WiUiams photon spectrum, the beamstrahlung spectrum 
and a backscattered laser spectrum. 
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Figure 20: Pseudorapidity cross sections dcr/drj^ measured in /ib per pseudorapidity unit, 
at the = 500 and 1000 GeV TESLA linear colliders [47]. Given are the distribu- 
tions for the Weizsacker-WiUiams photon spectrum, the beamstrahlung spectrum and a 
backscattered laser spectrum. 
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Figure 21: Transverse momentum cross sections da j dpi_ at the ^/s = 500 GeV TESLA 
hnear coUider [47]. Given are the distributions for the Weizsacker-WiUiams photon spec- 
trum, the beamstrahlung spectrum and a backscattered laser spectrum. Please note that 
the p_L distributions for the Weizsacker-WiUiams spectrum and the beamstrahlung spec- 
trum cross. At low p_L the beamstrahlung dominates, at high p_L the W.-W. spectrum 
dominates. 
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Figure 22: Decomposition of the jet cross section Jcrjet/c^-E'ijet as function of the jet 
transverse energy. The jets with -Eijet ^ 5 GeV were found from the PHOJET events using 
a cone algorithm. The jet acceptance region was restricted to |?yjet| < l-l- 
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Figure 23: Jet cross section Jcrjet/c^-E'ijet as function of the jet transverse energy. In order 
to study the influence of multiple soft and hard interactions in the model, we consider 
restricted cases with only single soft or single hard collisions. The curve with multiple 
soft, single hard collision coincides nearly with the full curve. The jets with -Eijet ^ 5 
GeV were found from the PHOJET events using a cone algorithm. The jet acceptance 
region was restricted to |?yjet| < 1.1. 
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Figure 24: Jet cross section Jcrjet/c^-E'ijet as function of the transverse energy. Here we give 
in addition the cross section with a visible energy cut (inside a calorimeter with < 2.1 
the cut is -^visible ^ 50 GeV). The jets with -Eijet ^ 5 GeV were found from the PHOJET 
events using a cone algorithm. The jet acceptance region was restricted to |?yjet| < 1.1. 
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Figure 25: Jet profiles in the pseudorapidity variable with and without the visible energy 
cut described in the caption of Fig. 24. 
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Figure 26: Cross section as function of the visible energy inside a calorimeter with < 
2.1. Only events with \/s^^ > 10 GeV are sampled from the improved Weizsacker- Williams 
spectrum. In a second curve we request jets with -Eijet ^ 5 GeV and impose a visible 
energy cut of 50 GeV. 



